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. Non-equilibrium fluid dynamics derived from the extended irreversible ther- 

I modynamics of the causal Miiller-Israel-Stewart theory of dissipative pro- 

Q I cesses in relativistic fluids based on Grad's moment method is applied to the 

' study of the dynamics of hot matter produced in ultra-relativistic heavy ion 

collisions. The temperature, energy density and entropy evolution are inves- 
tigated in the framework of the Bjorken boost-invariant scaling limit. The 
lO ! results of these second order theories are compared to those of first order 

I theories due to Eckart and to Landau and Lifshitz and those of zeroth or- 

Q^ ' der (perfect fluid) due to Euler. In the presence of dissipation perfect fluid 

dynamics is no longer valid in describing the evolution of the matter. First 
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• order theories fail in the early stages of evolution. Second order theories give 

a better description in good agreement with transport models. It is shown in 
which region the Navier-Stokes-Fourier laws (first order theories) are a rea- 
^ ' sonable limiting case of the more general extended thermodynamics (second 

order theories). 
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I. INTRODUCTION 



The study of space-time evolution and non-equilibrium properties of matter produced in 
high energy heavy ion collisions, such as those at the Relativistic Heavy Ion Collider (RHIC) 
at Brookhaven National Laboratory, USA and the Large Hadron Collider (LHC) at CERN, 
Geneva using relativistic dissipative fluid dynamics are of importance in understanding the 

observablcs. 

High energy heavy ion collisions offer the opportunity to study the properties of hot and 
dense matter. To do so we must follow its space-time evolution, which is affected not only 
by the equation of state but also by dissipative, non-equilibrium processes. Thus we need 
to know the transport coefficients such as viscosities, conductivities, and diffusivities. We 
also need to know the relaxation times for various dissipative processes under considera- 
tion. Knowledge of the various time and length scales is of central importance to help us 
decide whether to apply fluid dynamics (macroscopic) or kinetic theory (microscopic) or a 
combination of the two. The use of fluid dynamics as one of the approaches in modelling 
the dynamic evolution of nuclear collisions has been successful in describing many of the 
obscrvables [3,4]. The assumptions and approximations of the fluid dynamical models are 
another source for uncertainties in predicting the observables. So far most work has focused 
on the ideal or perfect fluid and/or multi- fluid dynamics. In this work we apply the rela- 
tivistic dissipative fluid dynamical approach. It is known even from non-relativistic studies 
[5] that dissipation might affect the obscrvables. 

The first theories of relativistic dissipative fiuid dynamics are due to Eckart [6] and to 
Landau and Lifshitz [7] . The difference in formal appearance stems from different choices for 
the definition of the hydrodynamical 4-velocity. These conventional theories of dissipative 
fluid dynamics are based on the assumption that the entropy four-current contains terms up 
to linear order in dissipative quantities and hence they are referred to as first order theories 
of dissipative fluids. The resulting equations for the dissipative fluxes are linearly related to 
the thermodynamic forces, and the resulting equations of motion arc parabolic in structure, 
from which we get the Fourier-Navier-Stokcs equations. They have the undesirable feature 
that causality may not be satisfled. That is, they may propagate viscous and thermal signals 
with speeds exceeding that of light. 

Extended theories of dissipative fluids due to Grad [8], Miiller [9], and Israel and Stewart 
[10] were introduced to remedy some of these undesirable features. These causal theories 
are based on the assumption that the entropy four-current should include terms quadratic 
in the dissipative fluxes and hence they are referred to as second order theories of dissipative 
fluids. The resulting equations for the dissipative fluxes are hyperbolic and they lead to 
causal propagation of signals [10,11]. In second order theories the space of thermodynamic 
quantities is expanded to include the dissipative quantities for the particular system under 
consideration. These dissipative quantities are treated as thermodynamic variables in their 
own right. 

A qualitative study of relativistic dissipative fluids for applications to relativistic heavy 
ions collisions has been done using these first order theories [12-17]. The application of sec- 
ond order theories to nuclear collisions has just begun [18-20], and the results of relativistic 
fiuid dynamics can also be compared to the prediction of spontaneous symmetry breaking 
results [21]. 



2 



The rest of the article is outUned as follows. In section 11 we review the basics of 
equilibrium thermodynamics and perfects fluids. In 111 the basic formulation of relativistic 
dissipative fluid dynamics will briefly introduced. In section IV we introduce the special 
hydrodynamical model of Bjorken. In section V we discuss the results of perfect fluids. In 
section VI we discuss the results of dissipative fluids. We flrst discuss the results of flrst 
order theories and then for the second order theories. In section Vll we summarize the 
results and discuss the need for hyperbolic theories for relativistic dissipative fluids. 

Throughout this article we adopt the units h = c = = 1- The sign convention used 
follows the time-hke convention with the signature (+,—,—,—), and if u°' is a time-like 
vector, u°'Ua > 0. The metric tensor is always taken to be g^" — diag(+l, — 1, — 1, — 1), 
the Minkowski tensor. Upper greek indices are contravariant, lower greek indices covariant. 
The greek indices used in 4-vectors go from to 3 {t, x, y, z) and the roman indices used 
in 3-vectors go from 1 to 3 (x, y, z). The scalar product of two 4-vectors , is denoted 
by Qfj^iy b" = a^^ bfj, = a • b. The scalar product of two 3-vectors is denoted by bold type, 
namely, a, b, a ■ b. The notations A^"^) = (^"^ + A^")/2 and A^""^^ = {A""^ - A^")/2 
denote symmetrization and antisymmetrization respectively. The 4-derivative is denoted by 
da = d/dx°'. Contravariant components of a tensor are found from covariant components 
by S-a/J^" = Af): 9ixa9ufsF°'^ = Ff,u and so on. 



II. BASICS OF EQUILIBRIUM RELATIVISTIC FLUID DYNAMICS 

The basic ingredients of equilibrium thermodynamics are the energy density e, the num- 
ber density n, and the entropy density s. The pressure p can be related to the other state 
variables. In this work these basic ingredients will be referred to as the primary thermody- 
namic variables. From the equation of state s = s{e,n), the temperature T and chemical 
potential emerge as partial derivatives: 

e+p^Ts + fxn , (1) 
Tds = d£ - ndn . (2) 

The chemical potential /i is the energy required to inject one extra particle reversibly into 
a flxed volume. The fundamental relation (1) then yields the last unknown thermodynamic 
function p. Next we introduce the thermal potential a and inverse-temperature /3 by the 
deflnitions 

, P^f ■ (3) 

Then equations (1) and (2) take the form 

s = (3{e + p) — an , (4) 
ds = Pde — adn . (5) 

In relativistic fluid dynamics it is advantageous to cast the thermodynamic quantities in 
terms of covariant objects, namely, the energy-momentum tensor T^", the 4-vector number 
current iV^ representing the net charge current, and the entropy 4-current S'^. For heavy ion 
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collisions, the conserved charge is usually taken to be the net baryon number. The covariant 
quantities are expressible in terms of a unique hydro dynamical 4-velocity u'^ (normalized as 
Uf^u^ = 1) and the primary state variables {e , n , p , s). 



N^^nu" , (6) 

T"" ^ eu'^u'' - pA''" , (7) 

S^'^su^' . (8) 

Here A^'^ = g^^ — u^u'^ is the spatial projection tensor orthogonal to m^, = 7(1, v) 



(7 = l/Vl — v^) is the (local) 4-velocity of the fluid, and gi^^ the Minkowski metric. The 
pressure can be related to the energy density and the number density by an equation of 
state. The equilibrium state will therefore be described by 5 parameters (e , n , u^). 

The thermodynamical relations can be written in a form that involves the covariant 
objects 5^, A^^ and T^'^ directly as basic variables: 

dSi" = -adNi" + (5AT>''' , (9) 

where 

= = ^ , (10) 

defines the inverse-temperature 4-vector. From (1) we have immediately 

S^" ^ pf3^' - aN" + P^T"'' . (11) 

It follows from (9) and (11) that 

dipP") = N^'da - T'^'dp^ . (12) 

Thus the basic variables A^'^, T^^'^ and can all be generated from partial derivatives of 
the fugacity 4-vector p(a, P)P^, once the equation of state is known. 

In the covariant formulation the hydrodynamical velocity is co-opted in a natural way as 
an extra thermodynamical variable. The covariant equations, expressed directly in terms of 
the basic conserved variables iV^ and T'^'^, can be manipulated more easily and transparently. 
Formulations for deviation from equilibrium can be carried out without much difficulty by 
making an assumption or requirement that (11) remains valid at least to first order in 
deviations (see the next section). 

The four-divergence of the energy-momentum tensor, the number 4-current, and the 
entropy 4-current vanish. That is, the energy and momentum, the number density, and the 
entropy density are all conserved. The covariant forms of particle, energy, momentum and 
entropy conservation are 



daN'' = ^ M"(9„n + n(9„M" = , (13) 

u^dfiT"^ = ^ u^daS + {e + p)do,u" = , (14) 

V^^aT'^^ = ^ (e + p)u^dpv'' - A'^^d^p = , (15) 

= ^ u^'das + sdau'^ = . (16) 
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Eq. (16) means that the entropy is conserved in ideal fluid dynamics. Note that s is constant 
along fluid world-lines, and not throughout the fluid in general. If s is the same constant on 
each world-line, that is, if Va-s = as well as s = 0, so that d^s = 0, then the fluid is called 
isentropic. Thus perfect fluids in equilibrium generate no entropy and no frictional type 
heating because their dynamics is reversible and without dissipation. For many processes in 
nuclear collisions a perfect fluid model is adequate. However, real fluids behave irreversibly, 
and some processes in heavy ion reactions may not be understood except as dissipative 
processes, requiring a relativistic theory of dissipative fluids. An equilibrium state is char- 
acterized by the absence of viscous stresses, heat flow and diffusion and maximum entropy 
principle, while a non-equilibrium state is characterized by the principle of nondecreasing 
entropy which arises due to the presence of dissipative fluxes. 

Perfect fluid dynamics has been successful in describing most of the observables [3,4,22]. 
The current status of ideal hydrodynamics in describing observables can be found in the 
latest Quark Matter Proceedings [1] and in [23,24]. Already at the level of ideal fluid 
approximation constructing numerical solution scheme to the equations is not an easy task. 
This is due to the nonlincarity of the system of conservation equations. Much work has been 
done in ideal hydrodynamics for heavy ion collision simulations (see e.g., [25]). 

III. NON-EQUILIBRIUM/DISSIPATIVE RELATIVISTIC FLUID DYNAMICS 

We give a covariant formulation of non-equilibrium thermodynamics. The central role 
of entropy is highlighted. Non-equilibrium effects are introduced by enlarging the space of 
basic independent variables through the introduction of non-equilibrium variables, such as 
dissipative fluxes appearing in the conservation equations. The next step is to flnd evolution 
equations for these extra variables. Whereas the evolution equations for the equilibrium 
variables are given by the usual conservation laws, no general criteria exist concerning the 
evolution equations of the dissipative fluxes, with the exception of the restriction imposed 
on them by the second law of thermodynamics. 

A natural way to obtain the evolution equations for the fluxes from a macroscopic basis 
is to generalize the equilibrium thermodynamic theories. That is, we assume the existence 
of a generalized entropy which depends on the dissipative fluxes and on the equilibrium 
variables as well. Restrictions on the form of the evolution equations are then imposed 
by the laws of thermodynamics. Prom the expression for generalized entropy one then 
can derive the generalized equations of state, which are of interest in the description of 
system under consideration. The phenomcnological formulation of the transport equations 
for the flrst order and second order theories is accomplished by combining the conservation 
of energy-momentum and particle number with the Gibbs equation. One then obtains an 
expression for the entropy 4-current, and its divergence leads to entropy production. Because 
of the enlargement of the space of variables the expressions for the energy-momentum tensor 
T^*^, particle 4-current A^^, entropy 4-current S'^, and the Gibbs equation contain extra 
terms. Transport equations for dissipative fluxes are obtained by imposing the second law 
of thermodynamics, that is, the principle of nondecreasing entropy. The kinetic approach is 
based on Grad's 14-moment method [8] . Por a review on generahzation of the 14-moment 
method to a mixture of several particle species see [26] . Por applications and discussions of 
the moment method in kinetic and transport theory of gases see e.g. [27] and for applications 
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in astrophysics and cosmology sec e.g. [28]. 

In the early stages of relativistic nuclear collisions we want to be able to describe phenom- 
ena at frequencies comparable to the inverse of the relaxation times of the fluxes. At such 
time scales, these fluxes must be included in the set of basic independent variables. In order 
to model dissipative processes we need non-equilibrium fluid dynamics or irreversible ther- 
modynamics. A satisfactory approach to irreversible thermodynamics is via non-equilibrium 
kinetic theory. In this work we will, however, follow a phenomenological approach. When- 
ever necessary we will point out how kinetic theory supports many of the results and their 
generalization. A complete discussion of irreversible thermodynamics is given in the mono- 
graphs [29-31], where most of the theory and applications are non-relativistic but includes 
relativistic thermodynamics. A relativistic, but more advanced, treatment may be found in 
[32-34] . In this work we will present a simple introduction to these features, leading up to a 
formulation of relativistic causal fluid dynamics that can be used for applications in nuclear 
collisions. 



A. Basic Features of Irreversible Thermodynamics and Imperfect Fluids 

The basic formulation of relativistic hydrodynamics can be found in the literature (see 
for example [7,35-37]). We consider a simple fluid and no electromagnetic fields. This fluid 
is characterized by, 

N'^(x) particle 4-current , (17) 
T^"{x) energy-momentum tensor , (18) 
S^{x) entropy 4-current , (19) 

where A = 1, ...,r for the r conserved net charge currents, such as electric charge, baryon 
number, and strangeness. N'^ and T^*^ represent conserved quantities: 

d,N^ = , (20) 
df.T^"' = . (21) 

The above equations are the local conservation of net charge and energy-momentum. They 
are the equations of motion of a relativistic fluid. There are 4 + r equations and 10 -|- 4r 
independent unknown functions. The second law of thermodynamics requires 

d^S'' > , (22) 

and it forms the basis for the development of the extended irreversible thermodynamics. 

We now perform a tensor decomposition of N'^, T'^", and with respect to an arbitrary, 
time-hke, 4-vector u^, normalized as m'^m^ = 1, and the projection onto the 3-space A'*'' 
orthogonal to u^. Note that 

A^"" = A"^ , A'^'u^ = , A/^" = A'^% A^ = 3 . (23) 

The tensor decomposition reads: 
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N^^nu^ + V^ , (24) 
T^'' = e u^v!" - p A^"" +2 W^" v!"^ + t^"" , (25) 
S^'^su^' + ^^ , (26) 



where we have defined 



W'' = q'' + hV'' , (27) 
^ T^M^ _ n A'^'' . (28) 

Here h is the enthalpy per particle defined by 

h^^-^. (29) 
n 

The dissipative fluxes satisfy the following orthogonality relations: 

Ui.V = 0, u^q'' = 0, Ui^W = 0, ^/^t'^'^ = 0, < = . (30) 

In the local rest frame (LRF) defined by -u^ = (1, 0) the quantities appearing in the 
decomposed tensors have the following meanings: 

(31) 
(32) 
(33) 

(34) 

(35) 
(36) 
(37) 
(38) 
(39) 

The angular bracket notation, representing the symmetrized spatial and traceless part of 
the tensor, is defined by 



n 






net density of charge 








net fiow of charge , 


e 






energy density , 


P + n 




— ''"A T^^'^ 
3 '^'^ 


pressure , 








energy fiow , 






W - 


heat fiow , 








stress tensor , 


S 






entropy density , 






A'^S'' 


entropy fiux . 



\ ( + A^^ A^) - i A>^^ A, 



(40) 



The space-time derivative decomposes into 

9'^ = ^^'^D + v^ «'^Vm = , (41) 

where 

D = u'^d^ convective time derivative , (42) 

yM = A'^^^i, gradient operator . (43) 
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In the LRF the two operators attain the given meanings. In this rest frame the projector 
becomes 



A^, = A^'^ = diag(0, -1, -1, -1), At = diag(0, 1, 1, 1) . (44) 

In the LRF the heat flow has spatial components only and the pressure tensor defined by 
pnv = ^ixrpar/Y^ j^jl) + tt'*'^ is also purely spatial. 

So far, is arbitrary. It has the following properties. Differentiating u^u'^ — 1 with 

respect to space-time co-ordinates, = {d/dx'^) = {dt, V), yields 

u^.d.u'' = , (45) 

which is a useful relation. 

There are two choices for [32]. It can be taken parallel to the particle fiux iV*, 



u 



I , (46) 

N^^nu'' . (47) 

This is known as the Eckart or particle frame and in this frame — 0. It can also be taken 
to be parallel to the energy flow. 



M^T'^'^ = eu^" . (49) 

This is known as the Landau and Lifshitz or energy frame and in this frame — 0. This 
implies that = —hV^. 

The two choices of velocity 4-flow have different computational advantage of each of the 
formulations. The Landau-Lifshitz formalism is convenient to employ since it reduces the 
energy-momentum tensor to a simpler form. The price for this is the implicit deflnition of 
the 4-velocity. The Eckart formalism has the advantage when one wants to have simple 
integration of particle conservation law. This choice is also more intuitive than that of 
Landau-Lifshitz. For a system with no net charge, the 4-velocity in the Eckart formalism is 
not well-deflned and therefore in general under this situation one should use the Landau- 
Lifshitz formalism. The Landau-Lifshitz formalism is also advantageous in the case of 
mixtures. 

Let us consider the Eckart deflnition of 4-velocity. With respect to this frame the particle 
number 4-current and energy-momentum tensor thus decomposes as follows: 

m = nu^' , (50) 
T^"^ = + 2g('^M^) + P^"" , (51) 

Pt^-^ = _(p + n) A'^" , A'^" TT^. = . (53) 

Successive terms on the right hand side of (53) represent the thermodynamical pressure, 
bulk stress and shear stress measured by an observer in the Eckart frame. 
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B. Conservation Laws 



We will now consider one type of charge, namely, the net baryon number. We also 
consider the Eckart frame for the 4-velocity definition. We insert the expressions for the 
number 4-current and the energy-momentum tensor in the conservation laws and project 
them onto the four velocity and the projection tensor. Using the orthogonality properties 
given in the previous section we obtain the following conservation laws. The equation of 
continuity (net charge conservation) d^N^ = 0, equation of motion (momentum conserva- 
tion) A'^dxT^^ = 0, and the equation of energy (energy conservation) u^d^T^^'^ = are, 
respectively, 

Dn = -n u'' > (54) 

{e + p + n)D = x/f^ip + n) - a;; tt"" + tt'^^d 

-[A:;Dg'' + 2g(^V.^^'^^] , (55) 

-y^q^ + 2q^'Du^ . (56) 

There are 14 unknowns, n, e, IT, q^, n^" , and u'^ and 5 equations. To close the system of 
equations we need to obtain 9 additional equations (for dissipativc fluxes) in addition to the 
5 conservation equations (for primary variables) we already know. The derivation of the 9 
equations will be presented in the following sections. 



C. Entropy 4-current and the generalized Gibbs equation 

Prom the phenomenological treatment of deriving the 9 additional equations we need 
the expression for the out-of-equilibrium entropy 4-current. The off-equilibrium entropy 
4-current S^'{N^',T^''') takes the form [10] 

= p{a, 13)13^ - aN" + P^T"" + Q^iSN", ST"", ...) (57) 

where Q'^ is a function of deviations SN^^ and ST'^" from local equilibrium 

ST^"" = T^"" - , SN^" ^N^"- , (58) 

containing all the information about viscous stresses and heat flux in the off-equilibrium 
state. 

Since the equilibrium pressure is only known as a function of the equilibrium energy den- 
sity and equilibrium net charge density we need to match/fix the equilibrium pressure to the 
actual state. We do this by requiring that the equilibrium energy density and the equilib- 
rium net charge density be equal to the off-equilibrium energy density and off-equilibrium 
net charge density. In Eckart 's formalism this is equivalent to 

ST^^'ui.u^ = SN^'ui,. (59) 

With the help of the expression for the divergence of pP" 
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d^ipPn = KA^ - Te79,Pu , (60) 

and the conservation laws for A'"'* and for T^^'^ the generalized second law of thermodynamics 
becomes 

d.S^' = -{6N>^)d,a + STi^^d^P, + d^Q'' . (61) 

Once a detailed form of Q'^ is specified, hnear relations between irreversible fluxes {6N'^, 
ST'^'^) and gradients {d(^^j3y), d^a) follow by imposing the second law of thermodynamics, 
namely, that the entropy production be positive. The key to a complete phenomenological 
theory thus lies in the specification of Q^. 



D. Standard Relativistic Dissipative Fluid Dynamics 

The standard Landau-Lifshitz and Eckart theories make the simplest possible assumption 
about Q^. That it is linear in the dissipative quantities (11, q^^i:^^). In kinetic theory this 
amounts to Taylor expanding the entropy 4-current expression up to first order in deviations 
from equilibrium. This leads to an expression of entropy 4-current which is just a linear 
function of the heat flux. 

This can be understood as follows: Take a simple fluid with particle current A^^. Let us 
choose (3^ = u'^/T parallel to the current of the given off-equilibrium state, so we are in 
the Eckart frame. Projecting (57) onto the 3-space orthogonal to gives 

= a;;^'^ = pq>' + a;; , q"^ m^t^^a^ , (62) 

so that 

^'^ = (f/T + second order terms , (63) 

which, to linear order, is just the standard relation between entropy flux $^ and heat flux 
q^. From (62) this implies that the entropy flux is a strictly linear function of heat flux 
q^, and depends on no other variables; also that the off-equilibrium entropy density depends 
only on the densities e and n and is given precisely by the equation of state s = Scq{e, n). 

Alternatively we may begin with the ansatz for the entropy 4-current S^: In the limit 
of vanishing 11, g^, and tt^'^ the entropy 4-current must reduce to the one of ideal fluid. The 
only non-vanishing 4-vector which can be formed from the available tensors m^, q'^ and tt'^'^ 
is Pq'^ where /3 is arbitrary but it turns out to be nothing else but the inverse temperature. 
Thus the first order expression for the entropy 4-current in the Eckart frame is given by 

S^ = su'' + ^ , (64) 

and one immediately realizes that 

^'^Y ' ^^^^ 

is the entropy fiux. Using the expressions for A'"'*, T^'^ and in the second law of thermo- 
dynamics d^S^ > and using the conservation laws Ui^d^/T^^" — 0, d^N^^ — and the Gibbs 
equation 
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the divergence of (64) gives the following expression for entropy production 



(66) 



(67) 



Notice that the equilibrium conditions (i.e., the bulk free and shear free of the flow and 
the constancy of the thermal potential,i.e., no heat flow) lead to the vanishing of each 
factor multiplying the dissipative terms on the right, and therefore lead to daS'^ — 0. The 
expression (67) splits into three independent, irreducible pieces: 



where 



X 
Xi" 



T 



-Du>' , 



(A^ A"; + A^f) - - A"" A, 



(68) 

(69) 
(70) 

(71) 



The angular brackets extract the purely spatial, trace-free part of the tensor. From the 
second law of thermodynamics, d^S^^ > 0, we see that the simplest way to satisfy the bilinear 
expression (67) is to impose the following linear relationships between the thermodynamic 
fluxes n, Q''', n^^'^ and the corresponding thermodynamic forces XJfjU'^, In T — Du'^, \j^'^u^'> : 



XT 



-Du" 



XnT^ 
e + p 



(72) 
(73) 
(74) 



That is, we assume that the dissipative fluxes are linear and purely local functions of the 
gradients. We then obtain uniquely, if the equilibrium state is isotropic (Curie's principle), 
the above hnear expressions. 

These are the constitutive equations for dissipative fluxes in the standard Eckart the- 
ory of relativistic irreversible thermodynamics. They are relativistic generalizations of the 
corresponding Newtonian laws. The linear laws allow us to identify the thermodynamic co- 
efficients, namely the bulk viscosity ({e,n), the thermal conductivity X{e,n) and the shear 
viscosity ri{6,n). 

Given the hnear constitutive equations (72) - (74), the entropy production rate (67) 
becomes 



2rjT 



QM 



+ 



> 



which is guaranteed to be non-negative provided that 

C>0, A>0, ?7>0. 



(75) 



(76) 
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Note that q^q^j, < which can be most easily proven from q^u^ = in the LRF. 

Using the fundamental thermodynamic equation of Gibbs the entropy evolution equation 
can be written in the following convenient form: 



(77) 



which can be found with the help of the fluid kinematic identity: 
where the fluid kinematic variables are given by 



(78) 




4-acceleration of the fluid , 



vorticity tensor , 
expansion tensor , 
volume expansion 



shear tensor . 



(79) 



The quantities defined here are the fluid kinematic variables. 

The Navier-Stokes-Fourier equations comprise a set of 9 equations. Together with the 5 
conservation laws d^T^'^ = d^N^ = 0, they should suffice, on the basis of naive counting, to 
determine the evolution of the 14 variables T'^'^ and N'^ from suitable initial data. Unfortu- 
nately, this system of equations is of mixed parabolic-hyperbolic-elliptic type. Just like the 
non-relativistic Fourier-Navier-Stokes theory, they predicts infinite propagation speeds for 
thermal and viscous disturbances. Already at the non-relativistic level, the parabolic char- 
acter of the equations has been a source of concern [38]. One would expect signal velocities 
to be bounded by the mean molecular speed. However in the non-relativistic case wave- 
front velocities can be infinite such as the case in the tail of Maxwell's distribution wich has 
arbitrarily large velocities. However, a relativistic theory which predicts infinite speeds for 
causal propagation contradicts the foundation or the basic principles of relativity and must 
be a cause of concern especially when one has to use the theory to explain observables from 
relativistic phenomena or experiments such as ultra-relativistic heavy ion experiments. The 
other problem is that these first order theories possess instabilities: equilibrium states are 
unstable under small perturbations [11]. 

Most of the applications of dissipative fluid dynamics in relativistic nuclear collisions 
have used the Eckart/Landau-Lifshitz theory. However, the algebraic nature of the Eckart 
constitutive equations leads to severe problems. Qualitatively, it can be seen from (72) 
- (74) that if a thermodynamic force is suddenly switched off/on, then the corresponding 
thermodynamic flux instantaneously vanishes/appears. This indicates that a signal propa- 
gates through the fluid at inflnite speed, violating relativistic causality. This is known as 
a paradox since in special relativity the speed of light is flnite and all maximum speeds 
should not be greater than this speed. This paradox was flrst addressed by Cattaneo [38] 
by introducing ad hoc relaxation terms in the phenomenological equations. The resulting 
equations conforms with causality. The only problem was that a sound theory was needed. 
It is from these arguments that the causal extended theory of Miiller, Israel and Stewart 
was developed. 
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E. Causal Relativistic Dissipative Fluid Dynamics 



Clearly the Eckart postulate (64) for Q'^ and hence S'^ is too simple. Kinetic theory 
indicates that in fact is second order in the dissipative fluxes. The Eckart assumption, 
by truncating at first order, removes the terms that are necessary to provide causality, 
hyperbolicity and stability. 

The second order kinetic theory formulation of the entropy 4-current, sec [8], was the 
starting point for good work during the 1950's on extending the domain of validity of con- 
ventional thermodynamics to shorter space-time scale. The turning point was Miiller's 
1967 paper [9] which, for the first time, expressed in terms of the off-equilibrium forces 
(n, q, vr*-' ) and thus linked phenomenology to the Grad expansion [8]. This marked the birth 
of what is now known as extended irreversible thermodynamics [29-31]. 

For small deviations, it will suffice to retain only the lowest-order, quadratic terms in 
the Taylor expansion of Q^^ , leading to hnear phenomenological laws. The most general 
algebraic form for that is at most second order in the dissipative fluxes gives [10] 



T 



T 



+ 



f 



(80) 



where /3^(£, n) > are thermodynamic coefficients for scalar, vector and tensor dissipative 
contributions to the entropy density, and aA{£, n) are thermodynamic viscous/ heat coupling 
coefficients. It follows from (80) that the effective entropy density measured by comoving 
observers is 



±,(^P,U'-P^q,q>^ + P27T,^.n^'') , 



(81) 



independent of aA- Note that the entropy density is a maximum in equilibrium. The 
condition m^Q^ < 0, which guarantees that entropy is maximized in equilibrium, requires 
that the Pa be nonnegative. The entropy flux is 



(82) 



and is independent of the Pa- 

The divergence of the extended current (80) together with the Gibbs equation (66) and 
the conservation equations (54), (55) and (56) lead to 



Tdf.s" = -n 



-q' 



(83) 



The simplest way to satisfy the second law of thermodynamics, 9^5''* > 0, is to impose, as in 
the standard theory, linear relationships between the thermodynamical fluxes and extended 
thermodynamic forces, leading to the following constitutive or transport equations. 
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7-nn + n = -co - 



\CT8, fe"" n 



-To V/. n - Ti Vi. tt); , 



rjTd, 



Here the relaxational times ta(£, n) are given by 



^2 „.A 



2r]T 



U TT, 



and the couphng coefficients Ti{s,n) are given by 

To = CcKo , Ti = XTai , T2 = 2r)ai . 



(84) 



(85) 



(86) 



(87) 



(88) 



A key quantity in these theories is the relaxation time r of the corresponding dissipative 
process. It is a positive-definite quantity by the requirement of hyperbohcity. It is the time 
taken by the corresponding dissipative fiux to relax to its steady state value. It is connected 
to the mean collision time tc of the particles responsible for the dissipative process, but the 
two are not the same. In principle they are different since r is a macroscopic time, although 
in some instances it may correspond just to a few tc- No general formula linking r and tc 
exists; their relationship depends in each case on the system under consideration. 

Besides the fact that parabolic theories are necessarily non-causal, it is obvious that 
whenever the time scale of the problem under consideration becomes of the order of or 
smaller than the relaxation time, the latter cannot be ignored. Neglecting the relaxation 
time in this situation amounts to disregarding the whole problem under consideration. 

Even in the steady-state regime the descriptions offered by parabolic and hyperbolic 
theories might not necessarily coincide. The differences between them in such a situation 
arise from the presence of r in terms that couple the vorticity to the heat flux and shear 
stresses. These may be large even in steady states where vorticity is important. There are 
also other acceleration coupling terms to bulk and shear stresses and the heat flux. The 
coefficients for these vanish in parabolic theories, but they could be large even in the steady 
state. The convective part of the time derivative, which is not negligible in the presence 
of large spatial gradients, and modifications in the equations of state due to the presence 
of dissipative fluxes also differentiate hyperbolic theories from parabolic ones. However, it 
is precisely before the establishment of the steady-state regime that both types of theories 
differ more signiflcantly. Therefore, if one wishes to study a dissipative process for times 
shorter than r, it is mandatory to resort to a hyperbolic theory which is a more accurate 
macroscopic approximation to the underlying kinetic description. 

Provided that the spatial gradients are not so large that the convective part of the time 
derivative becomes important, and that the fluxes and coupling terms remain safely small, 
then for times larger than r it is sensible to resort to a parabolic theory. However, even in 
these cases, it should be kept in mind that the way a system approaches equilibrium may be 
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very sensitive to the relaxation time. The future of the system at time scales much longer 
than the relaxation time once the steady state is reached, may also critically depend on r. 

The crucial difference between the standard Eckart and the extended Israel-Stewart 
transport equations is that the latter are differential evolution equations, while the former 
are algebraic relations. The evolution terms, with the relaxational time coefficients ta, are 
needed for causality, as well as for modelling high energy heavy ion collisions relaxation 
effects are important. The price paid for the improvements that the extended caiisal ther- 
modynamics brings is that new thermodynamic coefficients are introduced. However, as is 
the case with the coefficients (, A and rj that occur also in standard theory, these new coef- 
ficients may be evaluated or at least estimated via kinetic theory. The relaxation times ta 
involve complicated collision integrals. They are usually estimated as mean collision times, 
of the form r ~ ^— , where o" is a collision cross section and v the mean particle speed. 

The form of transport equations obtained here is justified by kinetic theory, which leads 
to the same form of the transport equations, but with extra terms and explicit expressions for 
transport, relaxation and coupling coefficients. With these transport equations, the entropy 
production rate has the same non-negative form (75) as in the standard theory. In addition 
to viscous/ heat couplings, kinetic theory shows that in general there will also be couplings 
of the heat flux and the anisotropic pressure to the vorticity. These couplings give rise to 
the following additions to the right hand sides of (85) and (86), respectively: 

+rqUJf,„q'' and + T^nf^uj„)x ■ (89) 

In the case of scaling solution assumption in nuclear collisions these additional terms do not 
contribute since they vanish. Also the resulting expression for d^S^ in general contains terms 
involving gradients of and [3^ multiplying second order quantities such as the bilinear 
terms {d^ai)qx'n'^'^ and {d^aQ)q^Il. In the present work where we will assume scaling solution 
these terms do no contribute to the overall analysis. 

It is also important to remember that the derivation of the causal transport equations is 
based on the assumption that the fluid is close to equilibrium. Thus the dissipative fluxes 
are small: 

|n| «p, {71^^.71^" «p, i-q^q'-f^ « e . (90) 

These conditions will also be useful in guiding us when we discuss the initial conditions for 
the dissipative fluxes. Considering the evolution of entropy in the Israel-Stewart theory, 
equation (77) still holds. 

F. Hyperbolicity, Causality and Characteristic Speeds 

Let us write the system of 14 equations in full, that is, the 5 conservation equations and 
the 9 evolution equations for dissipative quantities, in one single system. We will write down 
the linearized system; that is, we assume that the gradients of the macroscopic variables may 
be treated as small quantities of the order of the deviations from equilibrium. Consequently, 
terms hke V^Du^ etc., can be omitted because they are then of second order in deviations 
from equilibrium. The system of 14 equations, in a general frame, is then written as: 
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Dn + n\7^u'' + ^^V^ = , (91) 
(e + p)Du^' + v'^Ip + n) + DW^' - v.tt'^" = , (92) 
De + {e+p) Vm u'" + V^.W'' = , (93) 

PoDU - ao Vm g'' + + Vm (u^ + ^V^) = , (94) 

P.Dq'' + ao V'^ n + ai tt'^'^ + - - Du^^ = , (95) 

piDn'''' - ai q"^ + —n"" - y^'^w"^ - - V^'' ^ . (96) 

2r) n 

The system of 14 equations may be written as a quasilinear system of 14 equations in 
the form 

M%^{U'^)dJj'' = /^([/^) (A, = 1, 14) , (97) 

where M'^^iU'^) and /^{W^) can be taken to be components of 14 x 14 matrices and 14 
vectors. The right hand side contains all the collision terms, and the coefficients MJ'^(C/'^) 
are purely thermodynamical functions. 

Let S be a characteristic hypersurface for the system (97) and let 0(a;") = be the local 
equation for E. Then satisfies the characteristic equation 



det 



M|^(a«0) =0 . (98) 



0(a;") is a 3-dimensional space across which the variables are continuous but their first 
derivatives are allowed to present discontinuities [daU^] normal to the surface =^ [daU^] = 
U^{da4>)- The characteristic speeds are independent of the microscopic details such as cross 
sections. To solve the characteristic equation (98) we consider a coordinate system x" 
chosen in such a way that at any point in the fluid the system of reference is orthogonal and 
comoving. If is a function of and x^ only, the characteristic speeds can be determined 
from 

det(^;M^° - M^^) = , (99) 
where v is the characteristic speed defined by 

V = -do(f)/di(f) . (100) 

The 14-component vector U has been split into a scalar- longitudinal 6- vector Ul — 

[e,n,Il,u^,q^,TT^^]; two longitudinal-transverse 3- vectors (corresponding to the two trans- 
verse directions of polarization); Ulti = [u"^ , q"^ , tc^'^] and Ult2 = [u^ , q^ , tt^^] and purely 
transverse 2-vector Ut = [7r^^,7r^^ — tt^^]. Equation (99) for v accordingly splits into one 
6th-dcgrcc and two 3rd-degree equations. The purely transverse modes do not propagate. 
Hiscock and Lindblom [11] proved - as anticipated - that this general scheme, when applied 
to first order theories, always yields wave-front speeds v that are superluminal. 

We will be studying the dynamics of a pion fiuid in the hadronic regime and a quark- 
gluon plasma fiuid in the partonic regime. It is therefore important to check if these systems 



16 



conform with the principle of relativity under small perturbation of the equilibrium state. For 
a quark-gluon plasma we consider a gas of weakly-interacting massless quarks and gluons. 
We also consider such a system to have a vanishing baryon chemical potential {/ib = 0). 
This implies also that the net baryon charge is zero (n^ = 0). The equation of state is given 
by p = £/3. For massless particles or ultra-relativistic particles the bulk viscosity vanishes. 

In the absence of any conserved charge the convenient choice of the 4-velocity is the 
Landau-Lifshitz frame. In this case the characteristic equations for the wave-front speeds 
becomes very simple. For the longitudinal modes we get only the fast longitudinal mode 
(associated with the true acoustical wave). The absence of heat conduction has as a con- 
sequence the disappearance of the slow longitudinal propagation mode (associated with 
thermal dissipation wave) . The phase velocity of the fast longitudinal mode is given by 

vl = . (101) 

Using ~ 3/4 X 1/p (see [10] for the coefficients and Pi), we see that 

vl = l. (102) 

If we are considering only the shear viscosity we will get the above result. The wave-front 
speed (signal speed) for the transverse plane wave is given by 

Using ~ 3/4 x 1/p, this reduces to 

4=1. (104) 

For a pion fluid with vanishing chemical potential yUjr = we have, for the fast longitudinal 
mode, the following expression for the wave-front speed 

2 _ lPo + P2 + Pop2{e + p)dp/de 

and the wave-front speed of the transverse plane wave is given by 

4 = -7T^ 7 ■ (106) 

Note that for the pion system we are in the relativistic regime. Then the equation of state is 
taken to be that of a non-interacting gas of pions only. The bulk viscosity does not vanish. 
We show the dependence of the wave-front speeds and the adiabatic speed of sound on 
temperature in Fig. 1. 
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FIGURES 




0.1 1 10 100 

z = m/T 



FIG. 1. The transverse (vr), longitudinal {vl) and sound (vg) phase velocities in a pion gas, as 
a function of z = m/T. m is the mass of pion. 

For completeness and as an illustration we show here the results of a nucleon-nucleon 
system at a finite chemical potential. In this system all modes are present, including the 
thermal mode. The results are shown in Fig 2. 




z = m/T 

FIG. 2. The transverse (vt), longitudinal {vi, V2) and sound (vg) phase velocities in a nucleon 
^clS^ ctS ct function oi z = m/T. m is the mass of nucleon. 
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Thus in all the systems we will consider here (see next section) causahty is obeyed. 

IV. RELATIVISTIC NUCLEUS-NUCLEUS COLLISIONS 

In heavy ion colhsions there is much interest in the mid-rapidity region because the 
equation of state is simple for vanishing net baryon charge. Based on the observation that 
the rapidity distribution of the charged particle multiplicity dNch/dy is constant in the mid- 
rapidity region [39] , that is, it is invariant under Lorentz transformation in the mid-rapidity 
region, it is reasonable to assume that all other quantities like number density, energy density 
and dissipative fluxes also have this symmetry. Thus these quantities depend on the proper 
time, r, and 

r = t/7 = tVl - = Vf^ - . (107) 

The longitudinal component of the matter velocity is assumed to scale as = zjt which 
includes a boost-invariance for the central region. This is the Bjorken [39] scaling solution 
assumption for high energy nuclear collisions. In nucleus-nucleus collision experiments a 
central rapidity plateau is observed [1]. Thus one would expect that the measured quantities 
from the central region will be boost-invariant. By introducing the above scaling, the boost- 
invariance is automatically built into the model, simplifying the hydrodynamical equations. 
The effect of the transverse hydrodynamical flow is to remove the scale invariance at radial 
distances of the order of the nuclear radius. 

All particles are assumed to originate from one point in space-time, that is, along the 
constant r hypersurfaces. This symmetry is assumed to be achieved at the time of ther- 
malization after which the system evolution is governed by the equations of relativistic fluid 
dynamics for given initial and boundary conditions. 

Covariantly, the four-velocity can be expressed as 

= -(i,0,0,^) = — . (108) 

T T 

One deflnes the spatial rapidity y as 

.4ln(i±i) . (109) 

The spatial proper-time, r is invariant imder boosts along the 2;-axis. With the deflnitions 
of spatial rapidity and proper time one can find the inverse transformations 

t — Tcoshy , (110) 
z — Tsinhy . (HI) 

Then the four-velocity can be written as 

li" = (cosh y, 0,0, sinhy) . (112) 

The transformation of the derivatives is 
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dt 
d_ 



cosh y — sinh y 
— sinh y cosh y 



d_ 

ld_ 
\rdy J 



(113) 



We also note that using the transformation of derivatives and the definition of the 4- 
velocity we can write 



(114) 

(115) 



The recipe given in this section will be used to simplify the equations of relativistic fluid 
dynamics in the next sections. In this work we consider a (l+l)-dimensional scaling solution 
in which we have one nonvanishing spatial component of the 4- velocity in a (3+1) space- 
time. 



V. IDEAL FLUID DYNAMICS 

Assuming local thermal and chemical equilibrium, hydrodynamical equations can be 

derived from the Boltzmann equation. The perfect fluid approximation is the simplest 
approximation. The relativistic hydrodynamical equations follow from the conservation of 
the baryon number, energy, and momentum. From the conservation of energy-momentum, 
d^T^'^ — 0, multiphed by u'^, the relativistic Euler equation follows 

d^e +{e + p)d^,u^ - . (116) 

Considering the energy-momentum conservation, d/j^T^^ — 0, we change the variables 
t, z in the conservation laws for one-dimensional longitudinal motion in the absence of the 
conserved charges, 

dtT^^ + d,T'^ ^ , (117) 
dtT^^ + d,T'' ^ . (118) 

to the new variables r, y. We assume only a longitudinal boost-invariant expansion, u'^ = 
x^/t, (Bjorken scenario [39]), as might be the case at RHIC and LHC energies. Then the 
coupled system of partial differential equations reduces to 

The second equation (120) implies that there are no forces between fluid elements with 
different y. Thus in the central rapidity region thermodynamic quantities don't depend on 
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y. Prom equation (119), together with the first law of thermodynamics, one derives the 
conservation of the entropy. 



? = -- ■ (121) 

OT T 

which can be immediately integrated to give 

ST = SqTq = constant , (122) 

at constant y. Equation (122) shows that the entropy density decreases proportional to r 
independent of the equation of state of the fluid. 

For an ideal ultra-relativistic gas, such as the non-interacting quark-gluon plasma, e = 3p 
and the evolution of the energy density, entropy density and temperature, can be determined 
easily: 

To \ 



T(t) = 71) j , (123) 
s{t) = so- , (124) 

T 

sir) = So (^j . (125) 

Here Tq, eo, Sq and Tq are the initial time, energy density, entropy density and temperature, 
respectively. They are determined by the time at which local equilibrium has been achieved. 




FIG. 3. Proper time evolution for temperature, energy density and entropy density in the 
Bjorken model for nuclear collisions (longitudinal expansion only). 
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The evolution of the energy density, the entropy density and the temperature are shown 
as a function of proper time r in Fig. 3. For (l + l)-dimensional scahng solution the scalar 
volume expansion decreases as ^ ~ t~^. Due to additional mechanical work performed 
the energy density decreases as £ ~ r~^/^ and consequently the temperature decreases as 

T ~ T-V3. 

The results of a hydrodynamical model depend strongly on the choice of the initial 
conditions. Therefore a reliable determination of the initial conditions is crucial. The initial 
conditions can be taken in principle from transport calculations describing the approach 
to equilibrium, such as the Parton Cascade Model (PCM) commonly known as VNI [40] 
which treat the entire evolution of the parton gas from the first contact of the cold nuclei 
to hadronization or Heavy Ion Jet Interaction Generator (HIJING) [41]. For initializing a 
hadronic state one can use Ultrarelativistic Quantum Molecular Dynamics (UrQMD) [42]. 

Another frequently used relation between the initial temperature and the initial time is 
based on the uncertainty principle [43] . The formation time r of a particle with an average 
energy (E) is given by r{E) ~ 1. The average energy of a thermal parton is about 3T. Hence, 
we find tq — l/(3To). However, the formation time of a particle, the time required to reach 
its mass shell, is not necessarily identical to the thermalization time [43]. Consequently, 
determination of the initial conditions is far from being trivial. However, if data for hadron 
production are available, such as at SPS, they can be used to determine or at least constrain 
the initial conditions for a hydrodynamical calculation of observables such as the photon 
spectra [44]. 

The hydrodynamical model presented above for illustration is certainly oversimplified. 
Transverse expansion cannot be neglected, especially during the later stages of the fireball, 
significantly changing the observables at RHIC and LHC. In reality the expansion of the 
system will not be purely longitudinal; the system will also expand transversally. 

It might not be justified to use an ideal fiuid; dissipation might be important too. Hence, 
the Euler equation should be replaced by the Navier-Stokes [16,17] or even higher order 
dissipative equations [18-20]. 

Under the simplifying assumption of an ideal fiuid, the hydrodynamical equations can 
be solved numerically using an equation of state and the initial conditions, such as initial 
time and temperature, as input. The final results depend strongly on the input parameters 
as well as on other details of the model, as in the simple one-dimensional case. Also, it is 
important to adopt a realistic equation of state, in particular for the hadron gas. 



VI. NON-IDEAL/DISSIPATIVE FLUID DYNAMICS 

In dissipative fluid dynamics entropy is generated by dissipation. The dissipative quanti- 
ties, namely, H, g'* and n^^ are not set a priori to zero. They are specifled through additional 

equations. Since we will be working with a baryon-free system (n = 0), a convenient choice 
of the reference frame is the Landau and Lifshitz frame. The number current, energy- 
momentum tensor and the entropy 4-current in this frame are given by 

Nf" = nu^ + 1/^ , (126) 
T^'^ = eu^u'' - (p + H) N""" +7r'^'' , (127) 
S^^su^-^hV^ , (128) 
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Taking the divergence of the entropy 4-current with the help of conservation of net baryon 
number, d^N^ = 0, conservation of energy, u^d^T^^'^ — together with the relativistic 
Gibbs-Duhem relation in the form 



1 

n 



\7> 



gives the equation for the entropy flow: 



Td,, 



T 



/i 



(129) 



(130) 



where s is the entropy density. While in the case of a perfect fluid the entropy is conserved, 
here we see that the right hand side of (130) represents entropy production. In the Landau 
and Lifshitz frame the energy and momentum conservation laws are, respectively. 



(£ + p + Il)Du>' = V''(P + n) - A(; Va T^"" + T^'-'^Du, . 
The energy-momentum tensor in the Landau-Lifshitz LRF is given by 



(131) 
(132) 



-'-LRF 



fe \ 

(p + n-7r/2) 

(p + n-7r/2) 

VO (p + U + tt) J 



(133) 



This satisfy Tl^ = e — 3{p + H), = n'^^Uy = 0. To study the dynamics of the system 
it is necessary to obtain the energy-momentum tensor as seen by an observer at rest with 
respect to the Minkowskian coordinates. To this end it is necessary to apply a boost in the 
longitudinal direction. Using — (coshy, 0, 0, sinhy) we have 



fWcosh^y-V 

Vi 



V W cosh y sinh y 



W cosh y sinh y \ 

V± 
Wsinh^y + P / 



(134) 



with W — e + V the effective enthalpy density, V — p + H + n the effective pressure and 
V± = p + H — 7i/2 the effective transverse pressure. For perfect fluids tt = 11 = 0. It is clear 
that the effects of viscosity is to reduce pressure in the longitudinal direction and increase 
pressure in the transverse direction. In the first order theory 



r 
4 1 

TT = --Tj- 
3 T 



(135) 
(136) 



In the second order theory 11 and tt have to be determined from the second order transport 
equations. 

Already in the first order theory regime the effects of various dissipative processes are 
clear. Bulk viscosity corresponds to the isotropic gradient and its effect is to lower the 
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pressure. Heat conductivity is associated with energy flow from a temperature gradient. 
Shear viscosity, which is associated with anisotropic velocity gradients, lowers the pressure 
along the gradient. In the second order theory regime there are additional terms such as 
relaxation and coupling terms. At RHIC the local longitudinal pressure scales with the 
inverse proper time between collisions. Thus the pressure in the longitudinal direction is 
reduced. This in turn will lead to a reduction in the longitudinal work. This implies more 
transverse energy. Transverse expansion will be accelerated. 



A. Standard Dissipative Fluid Dynamics 

In the first order theories the dissipative fluxes are given by 



n = -C^ , 




" nT ' 




_e+p_ 







(137) 
(138) 

(139) 



The current is induced by heat conduction, = —q^/h. 

The (l + l)-dimensional scaling solution implies that the thermodynamic quantities de- 
pend on r only. Thus the scaling solution and the relations dr/dx^ = and df{T)/dx^ = 
Un{df{T)/dT) (where /(r) represent thermodynamic variables such as temperature, chemical 
potential and dissipative fluxes) lead to 



n = -C- , 

T 

Of) ^ • 1 2 

TT = — 77— smh y , 
3 T 

zz 4 1 ,2 

TT = — 77— cosh y . 
3 T 



(140) 
(141) 
(142) 

(143) 



Note that in the energy and entropy equations the quantity that appears is tt"" — tt^^, which 
is independent of y. Equation (141) implies that there is no heat conduction in the scaling 
solutions. This is independent of the fact that ti = 0, another condition that also makes 
vanish. The scahng solution also implies that 



(e + p)W = , v^'P = , A^;(9„(7r^" - A'^^H) = . 



(144) 



Before we continue, let us note that in (N+1) space-time 5'^ = N + 1^ = and N 
is the dimension of space. The symmetrized spatial and traceless part of the hydrodynamic 
velocity gradient is in general 



This is traceless since in {N + 1) space-time A^ — N 



(145) 
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The r dependence of thermodynamic quantities can be determined from the entropy 
and energy equations. We will only consider a (l+l)-dimensional scaling solution in (3+1) 
dimensions since wc only have shear viscosity. The energy and entropy conservation becomes, 
in first order theory, 



ds e + p 
dr T 

Similarly, the entropy equation becomes 



4:7] + C 1 



ds ^ s 
dr T 
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(146) 



(147) 



These two equations form the basis for the first order theory results in this work. In terms 
of the ratio of the nondissipative term to dissipative term we can write the above equations 
as 



de e + p 

dr T 
ds 



s + p 

Rt 

s 



s 

dr T Rt 



where the ratio i?, associated with the Reynolds number in [13,51], is defined by 



R- 



{e -\- p)t Tst 



C 



(148) 
(149) 

(150) 



In the case of perfect fluid, R~^ and hence the right hand side of Eq. (149) vanish. Thus the 
entropy density in the ideal fluid approximation is proportional to r~^, since st is constant. 

For this exploratory study a simple equation of state is used, namely that of a weakly 
interacting plasma of massless u ,d ,s quarks and gluons. The pressure is given by p = 
e/3 — aT^ with zero baryon chemical potential. That is, /j, — 0, e — 3p or s — AaT^ , 
r] — bT^ and ( — , a, b — constant. The solution of Eq. (149) becomes 



T 


'To' 


1/3 1 




. T . 






'To' 


"I 




. T . 





1 + 



1 + 



b 


1 - 


To 


2/3- 




GaroTo 




. T . 




) 



Rq ^ 


1 - 


To' 


2/3" 


1 


2 




. T . 







(151) 



where Tq and Rq are the initial values of the temperature and the Reynolds number at the 
initial proper time t — Tq. Here 



2 V 90 



(152) 



is a constant determined by the number of quark flavors and the number of gluon colors. 
In the case of massless particles the bulk pressure equation (137) docs not contribute since 
the bulk viscosity is negligible or vanishes [36] . The only relaxation coefficient we need is 
which, for massless particles, is given by — 3/(4p). The shear viscosity is given by [52] 
T] — bT^ where 
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is a constant determined by the number of quark flavors and the number of gluon colors. 
Here Nf is the number of quark flavors, taken to be 3, and ag is the strong flne structure 
constant, taken to be 0.4-0.5. 

A nonvanishing makes the coohng rate smaller, as expected. If Rq^ = 0, Eq. (151) 
gives 



T 



1/3 

(154) 



as is expected for an ideal fluid. The role of R can be examined by rewriting energy and 
entropy equations as 

As is seen in these equations, the energy density and entropy density increase if i? < 1 and 
decrease if i? > 1. Therefore there is a peak in s (and hence in T and in e) ii Rq < 1, and 
no peak if it!o > 1- -R = 1 is the critical value for Reynolds number. One notes that at this 
value thermodynamics quantities do not change with time. R itself increases monotonically 
towards a limiting value during the expansion stage as shown in Fig. 4. 
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FIG. 4. The r dependence of temperature T, entropy density s and the inverse of Reynolds 
number R"^ in the scahng solution with s and 77 are proportional to T^. 

In Fig. 4 the r dependence of T, s, and with two typical values of ag is compared 
with the perfect fluid case. For both values of c^s, T, s and have peaks, while for the 
perfect fluid they monotonically decrease and monotonically decreases. In both cases 
the initial conditions are those determined from the uncertainty principle. This way of flxing 
initial conditions has a strong restriction on the value of the initial Reynolds number. It is 
seen in Fig. 5 that Rq^ will always be greater than one. This is the condition that gives 
rise to the peaks in temperature, energy density and entropy density. In the same Fig. 5 
we also show the region for the condition Rq > 1 va. the T^TQ-as plane for arbitrary initial 
conditions (Tq, tq) at a given ccg. 

The condition for decrease in the thermodynamical variables is 



Rq>1 or TqTq > 



lb 
3 a 



where the initial Reynolds number can be calculated from 
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(157) 
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The above condition and the initial conditions tq ~ 1/ (3To) implies b < a. It also determines 
the region where one might apply the Navier-Stokes-Fourier laws. 

In (1+1) dimensions the (l+l)-dimensional scaling solution imphes that the effective 
viscosity coefficient is only the bulk viscosity (. This is also the case with the (3+1)- 
dimensional solution in (3+1) dimensions. This is due to the relation \j^u^ = l\.^^dpU^ 
resulting from scaling solution. The (l+l)-dimensional scaling solution in (3+1) dimensions 
implies that the effective viscosity coefficient is 4?7/3 + 
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FIG. 5. The dependence of Rq and the region which satisfies i?o > 1 is shown in the Us-ToTq 
plane; where the curve is the condition i?o = 1- 
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FIG. 6. The r dependence of temperature T and entropy density s for arbitrary initial condi- 
tions and for initial conditions from uncertainty principle. 

In Fig. 6 we show how initial conditions affect the overall space-time evolution. In the 
top panel we show the results for arbitrarily chosen initial conditions. For this set of initial 
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conditions dissipativc effects are small. On the other hand the bottom panel with the initial 
conditions from uncertainty principle indicates that dissipative effects are significant. 



B. Causal Fluid Dynamics 

We have seen in the previous section the importance of Reynolds' number in determining 
the effects of dissipation. One of the mathematical advantages of the previous section is the 
direct connection between the Reynolds number and the initial conditions (Tq, Tq). This is 
because the first order theory does not have well-defined initial conditions for the dissipative 
fiuxes, and the latter are related to the thermodynamic forces by linear algebraic expressions. 

In this section we determine the dissipative fiuxes from the transport equations. In the 
Landau-Lifshitz frame the transport equations are still given by (72)-(74) but with slightly 
different heat coupling coefficients in the bulk and shear viscous pressure equations. Under 
the scaling solution assumption those coupling terms do not contribute to the dynamics of 
the system. 

The (1 + l)-scaling solution in (3 + l)-dimensions yields the following structure of trans- 
port equations 



dU 

a7 



9^ 



1 1 



'PoT ' 



(159) 
(160) 




(161) 



In the last of the above equations A^'^ — for < < 1 and zero otherwise (because 
of only one non- vanishing spatial component of the 4-velocity). 

For the (1 + 1)- dimensional Bjorken similarity fluid flow in (3+1) dimensions the energy 
equation (131) becomes 



de e + p 
dr T 



1 



-$ + n- 



1 







(162) 



where $ = n^^ — n^^ is determined from the shear viscous tensor evolution equation (161) 




(163) 



The equation for IT will not be needed because C = 0. 

In this subsection we will distinguish the perfect fluid, flrst order, and second order 
theories by the quantity $: 
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$ = perfect fluid , (164) 

4 

^ — -ri/r first order theory , (165) 

d$ $ 1/1 Id //32\\ 2 11 , ^ , 



second order theory , 
The Reynolds number is generally given by 

^ = ^ or 5. (167) 

R goes to infinity in the ideal fiuid approximation. 

Let us first discuss the initial condition for $. For an ideal fluid $ vanishes since there 
are no dissipative fluxes. For the flrst order theories the initial condition for $ is not well- 
deflned and is given by the initial conditions (Tq, Tq). For the second order theories we have 
well-defined initial condition for $ since the dissipative fiuxes are found from their evolution 
equations. 

In deriving the transport equations it is assumed that the dissipative fluxes are small 
compared to the primary variables {e, n, p). For shear flux we require that 

K%.]'^' = yf*«p • (168) 

In terms of $ this condition can be written as 

^ « yf p ■ (169) 

Another important quantity in determining the importance of viscous effects is the ratio of 
the macroscopic time scale to the microscopic time scale 

Rr = , (170) 



which in our case is given by 



Rr^llTr. (171) 



In first order theories the question of how much a particular dissipative fiux is generated 
as a response to corresponding thermodynamic/kinematic forces in nuclear collisions is gov- 
erned by the primary initial conditions (Tq, tq). That is, one just reads off the value of $0 
from the the hnear algebraic expression for $. We have seen that for values of the Reynolds 
number less than one, the thermodynamic quantities increase with time. This might be 
signalling the instability of the solution. Alternatively this might imply that we are using 
the first order theories beyond their domain of validity . The primary initial conditions can 
in principle be extracted from experiments. These in turn will give us the value of $o- This 
value of $0 will eventually determine how the thermodynamic variables evolve with time. 
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This is clearly understood by looking at the ratio of the pressure term to viscous term, 
namely, R as already discussed in the previous section. 

In the second order theories the question of how much a particular dissipative flux is 
generated as a response to corresponding thermodynamic/kinematic forces in the early stages 
of nuclear reactions is not trivial but interesting. In order to find the viscous contribution 
to the time evolution of thermodynamic quantities wc need to solve the differential equation 
for $. Therefore one has to determine the initial conditions for $. Although we don't know 
the exact form of the initial value for $ we will discuss the limiting cases. The first and most 
important limiting case is based on the assumption made when deriving the second order 
theory transport equations, namely, that the dissipative fiuxes be small compared to the 
primary variables. For the shear viscous fiux this means that the shear viscous stress-tensor 
must be small compared to the pressure. The value of $ will always be less than p, hence 
the initial value $o will always be less than pq. This has an interesting consequence: the 
initial Reynolds number is always greater than one. Thus in second order theory under these 
conditions there will be no increase (and hence no peaks) of thermodynamic variables with 
increasing time. However in general the thermodynamic quantities will decrease with time 
for as long as the condition 



is satisfied, which in the present case implies that $0 ^ 4po- However, values of $ greater 
than the pressure p leads to unphysical negative effective enthalpy. Unlike in the first order 
theories where it is not always possible to address this problem of negative effective enthalpy, 
in the second order theories we are guided by the limitations which are embedded in the valid 
application regimes of the theories, namely, the condition {tt^^tThi^Y^^ <^ p. This condition 
guarantees that there effective enthalpy is always positive. 

There are other two ways of determining the initial conditions for $. The first one is by 
using the existing microscopic models such as VNI [40], (HIJING) [41] and (UrQMD) [42] to 
extract the various components of n^^'^ from T^^'^. Since we are dealing here with a partonic 
gas VNI seems to be a good choice for the present work. We will use the results from 
the improved version of VNI [46] to fit our calculation in order to extract the initial value 
for $. Another way of determining the initial value for $ is to extract the initial value of 
the Reynolds number experimentally. Two of the most experimentally accessible quantities 
are the multiplicity per unit rapidity dN/dy and the transverse energy per unit rapidity 
di^x/dy. A detailed study for the initial and boundary conditions for dissipative fluxes is 
needed to fully incorporate these fluxes into the dynamical equations for the thermodynamic 
quantities. 

In this section we focus mainly on the results of second order theories and compare them 
to the first order theories and the perfect fluid results. As it is the main focus of this work 
we will look at the dynamics of the quark-gluon plasma or partonic gas. 

The energy equation (162) and the viscous stress equation (163) can be written as 



^ <e+p , 



(172) 




T 



(173) 



12ar ' 



2aT$ 




9t 



(174) 
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For a perfect fluid and a first order theory the energy equation (173) can be solved ana- 
lytically to give (123) and (151) In the first order theory we do not have any relaxation 
coefficients. Then equation (163) gives $ = (4?7/3)/r. Eqs. (123), (151) and the numerical 
solution to the second order equations (173) and (174) will be used to study the proper time 
evolution of temperature. The other thermodynamic quantities, namely, energy density and 
entropy density, are related to the temperature by the equation of state. It is important 
to show the entropy results due to the importance of entropy in the theory of irreversible 
extended thermodynamics and due to the fact that entropy is related to multiplicity. 

We start by showing the dependence of the temperature evolution on the initial value of 
In Fig. 7 we show this dependence for two values of Ug. The dependence of first order 
results on ag is clear from the results presented before. The dependence of the second order 
theory on «s is different. A large value of ag modifies the power dependence for r > 1 fm/c 
to be more or less the same as the power dependence for r < 1 fm/c. Thus the effect of 
are important for r > 1 fm/c as can be seen from the top right panel of Fig. 7 where a large 
value of as is shown. 

In studying the dependence of the results on the initial conditions for $ we have also 
included some unphysical choices for illustrative purposes. For 0.1 < $o/Po < 1 the choice 
of $0 is important, but below $o = O.lpo the equation for $ gives the same contribution to 
the evolution of thermodynamic quantities. 




T [fm/c] T [fm/c] 
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FIG. 7. The time evolution of temperature for different choices of <I>o for given primary initial 
conditions (Tq, tq), and two different values of a^. The curves for different values of $o are from 
the second order theory. The bottom panel is for $q"'^ = $q**. 

It is also tempting to choose the initial conditions for the second order theory to match 
what the first order theory predicts to be the initial value of $. Note the order of the curves 
on the bottom panel of Fig. 7 . The second order theory predicts larger deviations than the 
first order theory. This should be exactly the same picture if both theories are synchronized 
in a regime where both are valid. Unfortunately, it is not trivial to make the reverse match 
of initial conditions. This situation will arise in natural way when both first order theory 
and second order theories are applied in the situation where they are both valid, as we will 
see later. 

Under physical initial conditions the second order theory gives a Reynolds number that 
is always greater than one. This can be seen from Fig. 8 where for illustrative purposes 
we also include curves for unphysical initial conditions for Note that $o = 4po is the 
maximum value before the solutions becomes unstable. This is a critical value that gives 
a Reynolds number Rq = 1. As expected the first order theory gives i? < 1 at the same 
time. Throughout this work, unless otherwise stated so, we use the primary initial conditions 
based on the uncertainty principle as already discussed and as we will briefly describe below. 
Under this prescription of primary initial conditions, which might be relevant for RHIC and 
LHC, the first order theories are not suitable in describing the dynamics of thermodynamic 
quantities. On the other hand, the second order theories are suitable in describing the 
physical process happening at earlier times. An advantage of the extended irreversible 
thermodynamics, or second order theories, are their ability to be applicable over a wide 
range of regimes. However, for a different choice of initial conditions both theories might 
yield similar results, as we shall see. 
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FIG. 8. Proper time evolution of the inverse Reynolds number for different values of $o for the 
given primary initial conditions. 

In what follows we will try to get close to the conditions that are reahzed in the laboratory. 
We will consider scenarios close to those at RHIC and LHC We will use the most common 
primary initial conditions as discussed below. But first, we have to estimate the initial value 
of $ for these two scenarios. We will use the recent results from VNI calculations for the 
proper time evolution [46]. We will make a fit to the data points and extract the initial 
value of 

Au+Au at = 200 AGeV Au+Au at = 200 AGeV 
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FIG. 9. The proper time evolution of energy density. The data points are results from VNI 
simulations and the curves are fits. 



Even though the motive is to extract an initial condition for there is something 
interesting in Fig. 9. In this figure a comparison between the perfect fluid approximation, 
the first order theory, and the second order theory is clear. The kinetic theory result, of 
course, differs significantly from the perfect fluid dynamics result. The first order theory 
obviously fails terribly. The essential point, however, is that the second order theory is in 
good agreement with the VNI results. Due to the prehminary nature of VNI results we 
cannot yet claim perfect agreement between the two approaches. However, the fact that 
both have similar power laws is striking. In the beginning it looks like and then later 
on r^^/^ for the VNI results. One expects that when the full three dimensional problem 
is studied within the fluid dynamical approach we might have even better agreement. The 
fitted value of $o is found to be about 0.2po which is, of course, a physical value. The 
value of Us used is about 0.5. For all RHIC results presented here we will use the expected 
primary initial conditions with $ = 0.2pQ and = 0.5. For the LHC scenario we will use 
the expected primary initial conditions with $o = 0.3po and as = 0.4. 

As a benchmark both the fluid dynamical and cascade models have been solved numeri- 
cally for same primary initial conditions and equation of state. This is done for consistency. 
It is apparent that hyperbolic models performs better than the parabolic ones, in agreement 
with VNIBMS simulations. Also for energy density there is a peak in the parabolic model 
which is absent in hyperbolic model. This spurious unphysical result highlights the differ- 
ence between the parabolic and hyperbolic model in region of large gradients. We remark 
that the initial state under consideration presents very steep velocity gradients. Therefore 
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this is an ideal benchmark for testing fluid dynamical models against transport models. 
Comparisons of Navier-Stokcs-Fourier results with transport models were made in [47] with 
NSF failing terribly for smaller cross sections. In that particular study the NSF also brought 
in the problem of negative effective pressure. The transport results however gave a much 
better description. What is important however is that the second order theory seems to do a 
better job even in this case. The latest results on this latter point to be published elsewhere 
are still under investigation and comparison to previous work on the effective pressure of a 
saturated Gluon plasma [48] is done. 

The initial temperatures expected are Tq^SOO MeV at RHIC and Tq^IOOO MeV at LHC. 
In Figs. 10 and 11 the initial time Tq is estimated by using the uncertainty principle [43]: 
To • (-E')o ~ 1 where {E)q ~ 3Tq for massless particles. This results in ro=0.13 fm/c at RHIC 
and ro=0.07 fm/c at LHC. The initial value used for $, must be specified independently 
for the second order theory. These values are legitimate since the second order theory is 
based on the assumption that the dissipative fluxes are small compared to the primary 
thermodynamic variables, namely p, n, and e. The effect of dissipation is more pronounced 
at the very early stages of heavy ion collisions when gradients of temperature, velocity, etc., 
are large. This can be seen by comparing Figs. 10, 11 with Fig. 12. At late times the effect 
of dissipation diminishes greatly. For comparison, in Fig. 12 we take a constant initial time 
To =1.0 fm/c, which is the characteristic hadronic time scale. From Figs. 10 and 11 Euler 
hydrodynamics predicts the fastest coohng. At hadronization times of order of 4 fm/c the 
perfect fluid dynamics solution decreases by a factor of about three. The first order theory 
fails badly even for this case where we have a very high initial parton density. The first order 
theory significantly under-predicts the work done during the expansion relative to the Miiller- 
Israel-Stcwart and Euler predictions. Thus the temperature decreases more slowly with the 
inclusion of dissipative effects. This would lead to greater yields of photons and dileptons. 
Also, the transverse energy and momentum would be reduced as the collective velocities 
are dissipated into heat. The system takes longer to cool down. This will delay freeze-out. 
The entropy, s = 4aT'\ is enhanced. This is important because entropy production can be 
related to the final multiplicity. 

A legitimate question to ask is: Do we really want to synchronize the initial conditions 
for both ideal fiuid, first order and second order theories? Given some initial conditions 
we want to investigate the importance of second order theories as compared to first order 
theories and perfect fiuids. That is, if one is given a set of well known initial conditions from 
experiment we want to see which of the theories best describe the dynamics of the system. 
Given an observable and a set of primary initial conditions we would like to see whether the 
microscopic cascade models, the ideal fluid, the flrst order theory or second order theory 
best describes the evolution of the system. 

Let us now analyze the differences between the second order and first order theories. 
The first thing we notice is that the Eckart-Landau-Lifshitz theory predicts that at early 
times the temperature will rise before falling off. This is more pronounced when we have 
small initial times. Naively one would expect that the system would cool monotonically as 
it expands, even in the case of dissipation where energy-momentum is conserved. On the 
other hand, it is seen that for large initial times and high temperatures the three theories 
have a similar time evolution. As can be seen from Fig. 12, all three cases start at the same 
point and then fall off with time. The difference stems from the fact that in the second order 
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theory the transport equations of the dissipative fliixes describe the evolution of these fluxes 
from an arbitrary initial state to an equilibrium state. The first order theory, though, is just 
related to the thermodynamic forces which, if switched off, do not demonstrate relaxation. 
Hence they are sometimes referred to as quasi-stationary theories. As can be seen from 
Fig. 12, it is before the establishment of an equihbrium state that the two theories differ 
significantly. 




X [fnVc] X [fnVc] 



FIG. 10. The proper time evolution of temperature and entropy for given primary initial con- 
ditions (To, To) for a RHIC scenario. 
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FIG. 12. The proper time evolution of temperature for RHIC scenario and LHC scenario. The 
initial time is arbitrarily chosen. The first order and second order theories overlaps. 



C. Hadronic Fluid Dynamics 

So far our focus has been on the quark-gluon plasma where the composition of the parton 
fluid enters the description through the form of the conservation laws and the equation 
of state. Now we study the dynamics of a pion fluid. Pions are the lightest hadrons. 
They are produced in abundance in ultra-relativistic collisions compared to heavier hadrons, 
particularly in the central region. It is therefore important to study their influence on the 
expansion. If pions are produced by hadronization of quark-gluon plasma, then dissipation 
encountered during their subsequent expansion may change the observables. The expansion 
in the central region conserves pion energy and momentum. Since pions carry baryon number 
zero, their total number is not conserved. Therefore, we expect the equilibrium number 
density of pions in a given volume to vary with temperature. 

The equation of state is approximated by that of a massless pion gas. Thus the pressure 
is given hj p = aT"^ with a = Qhir'^/dO where Qh = 3 is the number of degrees of freedom. 
The energy density and entropy density are given by e = 3 a and s = 4 a respec- 
tively. The bulk viscous pressure equation does not contribute for massless particles, since 
C — > [36]. For the (1+1) -dimensional Bjorken-type hydrodynamics the heat term in the 
energy equation will not contribute. Thus we need only the shear viscous pressure for this 
presentation. The energy density evolution equation is determined by 

+ (176) 



where 




dT T T 



perfect fluid , (176) 
standard theory , (177) 

extended theory , 
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where rj = Tt^/{2(32). For massless particles [32 = 3/(4 p), and this is used in the expression 
for Tj^. The primary transport coefficients of a massless pion gas are not that well-known. 
For chiral pions the expressions for shear viscosity and thermal conductivity are given in 
[26]. We will estimate the shear viscosity from the mean coUision time of the pions. The 
mean time between collisions of pions moving at (t;) ~ 1 is given by 

= — , (179) 
an 

where n{T) = bT^ with b = 3C(3)/7r^ is the pion density and a ^ 1 fm^ is an effective cross 
section. The quantity r^T^ = l/(cr6) fm^ is roughly constant for temperatures T > 100 
MeV. The shear viscosity can therefore be represented by 

77 = /,r with fv = l^- ■ (180) 

3 b a 

Using the transport and thermodynamic properties outlined here the energy and trans- 
port equations can be written as 

dT IT 1 1 ^ 

dr 2\r Tdr) 9 r ^ ^ 

The energy equation can be solved analytically for the perfect fluid and the first order 
(provided rj is constant) cases. But since we want rj to depend on temperature or time one 
must then solve the equations numerically, or first find the temperature evolution as done in 
the previous section. In the case of the perfect fiuid and the first order theories the energy 
equation is readily solved. 



T{t) = To 
T{t) = To 
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perfect fiuid , (183) 
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first order . (184) 



Hydrodynamics applies to the intermediate stages of the evolution until the system even- 
tually freezes out of local equilibrium into freely streaming fragments as it expands. At 
freeze-out hydrodynamics must break down, since the pion collision rate drops below the 
local expansion rate. The system will be purely longitudinal as long as its transverse motion 
can be neglected. Once the transverse rarefaction wave reaches the center, however, the 
fiow becomes effectively three dimensional, so that the density diminishes as t^'^. Freeze 
out then occurs because the mean free path, which then increases as r^, grows faster than 
the expansion time scale \ n/n \— t/3. 

In Fig. 13 we show the r dependence of temperature for the three different perfect 
fluid, a flrst order theory of dissipative fluids and a second order theory of dissipative fluids. 
Here we assume that the pion gas is produced at hadronization of quark-gluon plasma at 
T — A fm/c. As expected, in this regime, with the given initial conditions, the first order 
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and second order theories converge. This convergence is faster with increasing cross sections. 
The effects of viscosity are small but non-negligible. 




T [fm/c] T [fm/c] 

FIG. 13. The proper time evolution of temperature for or o" = 1.0 fm2 and a = 3.0 fm^. The 
initial conditions (Tq, tq) are arbitrarily chosen. The first order and second order theories overlaps. 

In Fig. 14 we assume that the pion gas is formed at Tq ~ l/(3To). As we know by 
now, the difference between the three theories is noticeable and first order theories arc not 
suitable. We see here also that the convergence of first order theory results and second order 
theory results will occur for large cross sections. 
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FIG. 14. The proper time evolution of temperature for a = 1.0 fm^ and a = 3.0 fm^. The 
primary initial conditions are those from uncertainty principle. 

The presence of dissipation in heavy ion reactions will have profound effects on the 
space-time evolution of the system. The freeze-out will be delayed. Temperature and 
energy density decrease more slowly. Enhancement of entropy production will increase the 
production of observed particles since the two can be related. Since the system takes longer 
to cool this will lead to an enhancement in the production of thermal signals (dileptons and 
photons) . 
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VII. CONCLUSIONS 



In this work I have given a comprehensive exposition of the non-equihbrium properties 
of a new state of matter produced in heavy ion reactions. In doing so I presented some basic 
features of non-equilibrium fluid dynamics. I studied the space-time description of high 
energy nuclear collisions. I chose to follow the phenomenological approach and left kinetic 
theory approach for future work. 

The ultimate aim is to bridge the phenomenological theory with the kinetic theory of 
the matter produced in heavy ion collisions. In doing so I made use of the dissipative fluid 
dynamics. The connection between the macroscopic theory and microscopic theory enters 
through the transport coefficients of the matter. The equation of state provided closure to 
the system of conservation equations. 

I demonstrated that extended irreversible thermodynamics provides a consistent frame- 
work to simulate and study the space-time evolution of ultra-relativistic nuclear collisions, 
from some initial time to the flnal particle yield. Although this approach relies on a number 
of fundamental assumptions and is far from providing an accurate quantitative description, 
it has the advantage of wide applicability. 

The advent of accelerators such as RHIC and LHC provide an opportunity for studying 
the dynamics and properties of the matter at very high energy density. In the description 
of the evolution of such a system, it is mandatory to evaluate, as accurately as possible, the 
order of magnitude of different characteristic time scales, since their relationship with the 
time scale of observation will determine, along with the relevant equations, the evolution 
pattern. This is rather general when dealing with dissipative systems. It has been my 
purpose here, by means of simple model with simple equation of state and arguments related 
to a wide range of time scales, to emphasize the convenience of resorting to hyperbolic 
theories when dissipative processes, either outside the steady-state regime or when the 
observation time is of the order of or shorter than some characteristic time of the system, 
are under consideration. Furthermore, dissipative processes may affect the way in which the 
system tends to equilibrium, thereby affecting the future of the system even for time scales 
much larger than the relaxation time. 

In the early stages of heavy ion collisions, non-equilibrium effects play a dominant role. 
A complete description of the dynamics of heavy ion reactions needs to include the effects 
of dissipation through dissipative or non-equilibrium fluid dynamics. As is well-known, 
hyperbolic theories of fluid dissipation were formulated to get rid of some of the undesirable 
features of parabolic theories, such as acausality. It seems appropriate therefore to resort 
to hyperbolic theories instead of parabolic theories in describing the dynamics of heavy ion 
collisions. Thus in ultra-relativistic heavy ion collisions, where the fluid evolution occurs 
very rapidly, the second order theories, due to their universality, should be used to analyze 
collision dynamics. 

Unlike in flrst order theories, where the transport equations are just the algebraic rela- 
tions between the dissipative fluxes and the thermodynamic forces, second order theories 
describe the evolution of the dissipative fluxes from an arbitrary initial state to a flnal 
steady-state using the transport equations. The presence of relaxation terms in second or- 
der theories makes the structure of the resulting transport equations hyperbolic and they 
thus represent a well-posed initial value problem. 
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The consequences of non- ideal fluid dynamics, both first order (if apphcable) and second 
order were demonstrated here in a simple situation, that of scaling solution assumption and 
simple equation of state. A more careful study of the effects of the non-ideal fluid dynamics 
on the observables is therefore important. Conversely, measurements of the observables re- 
lated to thermodynamic quantities would allow us to determine the importance and strength 
of dissipative processes in heavy ion collisions. 

Here we have used only a simple equation of state in a simplified model of high energy 
nuclear collisions. A more realistic situation (including transverse expansion) will require 
careful analysis of both the transport coefficients and the equation of state which are em- 
ployed in the full set of the equations. It is then that one may have a better understanding 
of when to use either of these theories in the context of relativistic heavy ion collisions. 

A. The Need for Hyperbolic Theories of Relativistic Dissipative Fluid Dynamics for 

Nuclear Collisions 

The need for hyperbolic theories of dissipative fiuid dynamics has been emphasized before 
[53] . But it is in ultra-relativistic nuclear collisions where it needs to be stressed. 

The first order theories of dissipative fiuid dynamics have a long history and proved 
very useful, in some special cases in understanding observables. However, they exhibit 
some undesirable features, such as acausality, that prompted the formulation of hyperbolic 
theories of dissipative fiuid dynamics to get rid of them. This was achieved at the price of 
extending the set of field variables by including the dissipative fiuxes, such as the heat current 
and viscous stresses, on the same footing as the primary ones, such as energy densities, 
equilibrium pressures, and number densities, thereby giving rise to a set of more physically 
satisfactory but involved theories from the mathematical point of view. These theories have 
the additional advantage of being backed by a kinetic theory of gases (Grad's 14-moment 
approximation). Second order theories, being hyperbolic in structure, lead to well-posed 
initial-value (Cauchy) problems. They also lead to causal propagation. Unlike the first 
order theories, second order theories have relaxation terms which permit us to study the 
evolution of the dissipative fiuxes. 

In addition to causality violation, the Eckart theory leads to unstable equilibrium states 
under small perturbations. A dissipative hot and dense nuclear matter may very rapidly 
tend towards a quasi-stationary state that is adequately described by the Eckart theory. 
However, this needs to be checked a priori. The Eckart theory can describe steady-state 
regimes. However it is clearly unable to deal with the evolution towards these regimes, or 
with the overall dynamics of the fiuid, in a satisfactory way. 

It can be argued that the infinite speeds of propagation predicted by parabolic theories 
and the unphysical behaviour of the space-time evolution of thermodynamic quantities is 
just a consequence of using them outside their range of applicability [14]. However hyperbolic 
theories also share this constraint but they naturally predict speeds which are of the order 
of the characteristic speeds of the problem at hand. Therefore hyperbolic theories have a 
much wider range of applicability than parabolic theories. 

Aside from the relaxation time r, hyperbolic theories introduce a certain number of new 
coefficients that couple the different dissipative fiuxes. This is rather natural, since these 
theories are designed to explain more complex phenomena than parabolic theories do. These 
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quantities are not free parameters that one could choose at wilL These new coefficients are 
restricted by the convexity of the entropy function, and can be calculated explicitly with 
the help of kinetic theory or fluctuation theory. It turns out that they are functions of the 
equation of state. 

In summary, although parabolic theories have proved very useful for many practical 
purposes, they appear to fail hopelessly in describing the dynamics of heavy ion collisions. 
In contrast, hyperbolic theories successfully give a better description in agreement with 
transport models and hopefully they will be able to predict the experimental results. Thus 
hyperbolic theories are more reliable. In the steady-state, under the conditions mentioned 
before, and for times exceeding r both theories converge. Based on the results presented here 
one can only stress the convenience of using hyperbolic transport equations when parabolic 
theories either fail or the problem under consideration happens to lie outside the range of 
applicability of parabolic theories. 

I will conclude by pointing out some of the advantages and challenges of the second order 
theories. The challenge we face is the increase in the space of thermodynamic variables. We 
now have, in addition to the transport coefficients, new coefficients in the problem. These 
are the relaxation coefficients j3i and the coupling coefficients a^. These new coefficients 
depend on the primary thermodynamic variables, such as n, e and p and therefore are deter- 
mined by the equation of state. Like viscosity and thermal conductivity, which are required 
to be positive by the second law of thermodynamics, these new coefficients are constrained 
by hyperbolicity requirements. In principle , in order to solve the second order relativistic 
dissipative fluid dynamic problem, one still needs the equation of state, initial conditions 
and the transport coefficients. Another challenge involves solving the equations numerically. 
Finally one would like to compare the results of non-equilibrium fluid dynamics to observ- 
ables. This will require an effort to solve the full system of equations numerically. Efforts 
in this direction are under way [54] 

The study of non-ideal or non-equilibrium fluid dynamics will be important for construct- 
ing hydro- molecular dynamic schemes [49]. Hyperbolic theories might prove to be convenient 
in such scheme. A resulting hydro-molecular dynamic scheme can then be compared to the 
studies of multi-fluid dynamics [50] . 

VIII. OUTLOOK 

There are important questions that need to be investigated in order to tackle the chal- 
lenges faced by hyperbolic theories. An important question is the measurability of the 
dissipative fluxes. The heat flux through a system may be simply evaluated by measuring 
the amount of energy transported per unit area and time. The viscous pressure can be 
measured from the tangential shear force exerted per unit area. In practice, it may be dif- 
ficult to evaluate these quantities at each instant and at every point. From kinetic theory 
these fluxes can be simulated from microscopic transport models such as HIJING, VNI and 
UrQMD. 

A thorough study of transport coefficients is needed. The results from lattice QCD 
should give clear predictions for these coefficients. So far there are no reliable phenomeno- 
logical expressions for the coefficients. The shear viscosity results presented in [55] using 
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the microscopic model UrQMD serve as a starting point for future calculations of transport 
coefficients. 

A thorough study of the equation of state is required in order to be consistent with the 
non-equilibrium description of matter. In addition to the initial/boundary conditions for 
the primary /equilibrium variables a thorough study of initial/boundary conditions for the 
dissipative fluxes is necessary. This is required for the evolution equations of the fluxes. 
For example, one would like to know how much of a particular dissipative flux is generated 
as a response to an associated thermodynamics force in the early stages of the heavy ion 
collisions. This in turn gives us information about the initial entropy generated as a result 
of dissipation. 

For the complete description of the dynamics of viscous, heat conducting matter we need 
to consider more realistic situations: a system that expands in both the longitudinal and 
transverse directions [54] and we need a full (3+l)-dimensional solutions to the conservation 
and evolution equations. This will require extensive numerical computation. This is a 
challenging but interesting problem. In order to understand the observables we need a full 
formulation of hyperbolic theory that should be tested against other models/theories. 

To solve imperfect hydrodynamics amounts to knowledge of (i) a realistic equation of 
state, (ii) reliable transport coefficients (ii) realistic initial/boundary conditions and (iv) 
numerical computational algorithm. In order to understand the observables from RHIC and 
LHC knowledge of these requirements is needed and that is what I have started to do. I 
hope that this will gather momentum in due course. 
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